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The study gives a decisive answer to the recently risen question about the type and origin of 
interaction between spin vortices and spin vacancies in 2D spin models. The approach is based on 
the low-temperature approximation of the 2D XY model known as the Villain model and does not 
involve any additional approximations, thus preserving the lattice structure. The exact form of the 
Hamiltonian describing a system of topological charges and a vacant site supports the attractive 
type of interaction between the vacancy and the charges. 

The quantitative difference between the characteristics of the vortex behavior in the 2D XY and 
Villain models due to the different energy of the vortex "cores" in the two models is pointed out. 
This leads to a conclusion that the interaction between a vortex and a spin vacancy and between a 
vortex and the antivortex differs quantitatively for small separations in the two mentioned models. 

PACS numbers: 05.50.+q; 75.10 



I. INTRODUCTION 

The term spin vortex has become common in theoreti- 
cal and experimental studies of magnetic materials. It is, 
in fact, a particular case of a more general class of physi- 
cal/mathematical objects called topological defects [H, 0- 

Although, strictly speaking, topological defects can be 
defined only in terms of a continuous field, similar for- 
mations can be observed in lattice spin models. More- 
over, it is the spin vortices that are responsible for the 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
in the 2D XY model (or, speaking more generally, 
in classical 2D easy-plane magnets). 

Most of the theoretical studies of the vortex proper- 
ties are limited to the low-temperature continuum model 
proposed by Kosterlitz and Thouless [|| (KT model, here- 
after). However, this approach obviously cannot give 
satisfactory results, when essentially "discrete" phenom- 
ena, as the effects induced by a spinless site, are studied. 
The lack of theoretical studies regarding spin vortices 
on a lattice and the related problem of spin vortex-spin 
vacancy interaction is the principal motivation for the 
present work. 



A. Spin vortices 

The 2D XY model is usually defined as a system of 
two-component spins S r of unit length which states can 
be represented by a polar coordinate —tt < 8 < tt: S r = 
(cos^rjSin^,.), placed at sites r of a square lattice, and 



described by the Hamiltonian 

H 2 DXY = JY, [1 - COS(0 r - r ,)] 



(1) 



or 
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Close enough to the ground state we have 8 r — 6 
±2tt for neighboring sites r, r'. 

Generally, considering two neighboring spins at sites r 
and r', one can encounter the two situations: \8 r — T i \ < 
7T and \9 r — 8 r >\ > tt (the situation \6 T — 9 r >\ = tt can be 
neglected) . In order to define spin vortices in the system 
under consideration, let us introduce the lattice of sites 
R, dual to the original lattice (the dual lattice is the set of 
all the centers of elementary cells of the original lattice) , 
and consider only those bonds (R, R') which intersect the 
bonds (r, r') of the original lattice for which \9 r — 6 r '\ > tt. 

In order to consider the bonds of interest in a system- 
atic way, let us say that r = (x, y),r' = (x 1 , y') define the 
bond (r, r') if r' = (x + a,y) for a horizontal bond and 
r' = (x, y + a) for a vertical bond, where a is the lattice 
spacing. The same rule is imposed for bonds of the dual 
lattice. Now, we can ascribe to every bond (R, R') a di- 
rection defined by the sign of 9 r — 6 T > of the intersected 
bond (r,r'): (R — > R') if 8 r — 9 r > > tt for a horizontal 
bond (R, R') or 9 r i — 6 T > tt for a vertical bond (R, R'), 
and (R' — > R) in the opposite case. The introduced rep- 
resentation is unique for a given microstate of the spin 
system (the revers statement is not true, of course). 

The most basic structural unit that can be distin- 
guished in the representation we have built is a path L 
(cither straight or steps-like) connecting two sites of the 
dual lattice, formed by one or several bonds connected 
together so that their directions comply with some gen- 
eral direction of the path. In the most general case, that 
path can be either closed or not closed. 

While a closed path L represents a trivial situation, the 
spin configuration with L starting and ending at different 
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sites of the dual lattice is of great interest and is called a 
vortex-antivortex pair (it can be said that the vortex and 
the antivortex are centered at the ends of the path L). 
The above concerns vortices with topological charges ±1; 
pairs of vortices with higher values of topological charge 
can be defined in terms of several paths that start at the 
vortex and end at the antivortex. Paths that start at the 
same site but end at different sites correspond to clusters 
of vortices with different absolute values of charge (for 
example, +2 and —1, —1). 

The regions around the vortex origins are character- 
ized by significant disorientation of spins and are called 
"cores". Moderate spin- wave excitations, when 6 r — 8 T i ~ 
0,±27r everywhere except for the vortex "cores", cannot 
destroy the vortex-antivortex pair unless the two defects 
annihilate at the same point. 

Short-range exchange forces between spins lead to 
long-range effective interaction between vortices. The 
energy of this interaction can be explicitly singled out 
in the Hamiltonian of the Villain model, and turns out 
to depend only on the essentially inherent characteristic 
of the vortices called topological charge [fj. It will be 
shown that the logarithmic asymptotic form, obtained 
by Villain for the attraction energy of the vortex and the 
antivortex at large separations, in fact holds sufficiently 
well on a lattice up to the smallest possible separation of 
one lattice spacing a (if neglecting the subtle anisotropy 
effects). 

On the contrary, the corresponding energy which we 
estimate for the 2D XY model turns out to deviate from 
the logarithmic law at small separations. In particu- 
lar, our result for the energy needed to create a vortex- 
antivortex pair is approximately 6.6J (J is the coupling 
constant), in contrast to 9. 9 J of the Villain model @, 
and in reasonable agreement with the recent results of 
Monte Carlo simulations @, Q ■ 

The details of the results announced above can be 
found in Section HT1 



in the continuous spin field. Subsequently, the same au- 
thors denied this nonphysical result, on the basis of their 
spin dynamics simulations [Til ]. 

The problem was resolved phcnomcnologically, postu- 
lating that the vacancy does not change the vortex struc- 
ture (or the change is negligible) [ljj. Under this as- 
sumption, the KT theory led to the attractive interac- 
tion which agreed with the results of computer simula- 
tions. However, this approach, giving correct qualitative 
picture, was not able to describe the particular details of 
the lattice under consideration. 

In our study, based on the Villain model, we obtain 
the effective Hamiltonian describing interaction between 
spin vortices and spin vacancies on a square lattice. 

For example, as it will be shown in this paper, the 
interaction energy for an individual vortex of topological 
charge q at point R and a spin vacancy at r reads: 

S(|r-R|) = -( 7 r-l)^^ + 0(|r-R|- 2 ) , (2) 

i.e. the vacancy and the vortex attract each other. 

Eq. ([2]) is the asymptotic expression which in fact holds 
well enough for separations as small as just a few lattice 
spacings. It will be argued that in the 2D XY model 
this energy differs considerably from @ for small sepa- 
rations |r — R|. For example, the vortex-on- vacancy pin- 
ning energy of the Villain model E{a) = — (3tt — 4) Jq 2 ~ 
— 5.425J<7 2 , in contrast to that of the 2D XY model ob- 
served in spin dynamics simulations, —3. 54 J [TTj . and 
other numerical studies, —3. 178 J [Tol ]. 

The details of the results announced here can be found 
in Section [TTT1 



II. VORTICES IN THE VILLAIN AND 2D XY 
MODELS 

A. Topological charges in the Villain model 



B. Spin vacancies 

An aspect of the spin vortex behavior, which only re- 
cently drew attention of the researchers, is the effective 
interaction with nonmagnetic inclusions in the lattice 0- 
fTTj . Such spin vacancies are part of the models with 
quenched disorder [T2T4TH and the lattice gas spin mod- 
els [Tfllljj]. Here, we will focus, however, not on the ther- 
modynamic quantities, but on the effective Hamiltonian 
which describes the interaction between spin vortices and 
vacancies. 

To our knowledge, the first theoretical works devoted 
to this problem demonstrated global deformation of the 
vortex structure caused by a single vacancy and repul- 
sive interaction between the vortex origin and the va- 
cancy [T^ . This result was essentially caused by an 
application of the KT continuum model which required 
representation of the vacancy by a cutout of a finite size 



Studying the low-temperature properties of the model 
([T]) it would be natural to apply the spin- wave (harmonic) 
approximation (SWA), i.e. to replace 1 — cos(# r — 8 r i) in 
(UJ with ^(0 r —9 T ') 2 . Indeed, this allows to examine many 
important pro perties of the low-temperature phase of this 
model jl8T - l20l ]. However, the states with \9 r — 8 r '\ > 
7r, which are crucial when considering spin vortices, will 
have non-physical energy in this case. So, the proper 
harmonic approximation must be 

H 2 DXY -2^2 6r ' ~ 2lTm (0r ~ 8 r ,) f (3) 

(r,r'> 

with 

r +1, 9 T -0 T . > tt; 
m(6 r -6 T ') = { -1, r - 6 T , < -tt; 

{ o, \e r -e r ,\ < tt. 
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At low temperatures, m(9 r — 9 r >) can be considered as 
independent degrees of freedom taking discrete values 
0, ±1. In turn, this leads to the Hamiltonian of the Vil- 
lain model: 



»4e 



— 6 r i — 27TTO r 



(4) 



(obviously, m r y - 
Assuming that 



-m T ' tT ). 



(5) 



where (p r and ip r are chosen so that \tp r — (p r i\ < n for 
any pair of spins in the system, i.e one can say that the 
field if r is vortexless, and all the vortices are "contained" 
in tp r , the Hamiltonian (j4|) can be written as: 



H = — ^ [(<Pr - Vv') 2 + (ipr ~ i>r' ~ 2TTm r ^) 2 ] 



<r,r'> 



J*^2<Pr^(lpr - V'r+u - 2lTm r , T+u ) 



(6) 



with u = (±o, 0), (0, ±a) and lattice spacing a. 

Following Villain jf|, one can choose ipr{{m r ,T'}) such 
that ip r and ip r decouple in the Hamiltonian, i.e. the last 
term in ([6]) vanishes: 



(ip T — i^r+u — 27rm r . r + u ) = for all r 



(7) 



This is realized when 



^Y.{ ( m3 < 4 _ m h2)Isc(x -X,y-Y) (8) 

Ft, 

+ (m4,i - m 2 ,3)Isc(y — Y,x — X) + {mi, 2 
- ffi2,3 + rn 3A - m i ,i)I ss {x — X,y — Y) j 



(see Fig. [T]) where R are sites of the dual lattice, which 
are situated in the centers of elementary cells of the orig- 
inal lattice, and functions I sc and I ss are given by (|A1[) 
and (|A2[) . (The asymptotic properties of I sc and I ss are 
analyzed in Appendix [A"l) In fact, Eq. (|SJ) is another way 
of presenting the expression obtained by Villain Q . 

In © , the sum over R = {X, Y) spans the sites of the 
dual lattice, while coordinate r represents a site of the 
original lattice, therefore, X — x and Y — y can be always 
presented as {2n — 1)# , where n is an integer. The short 
notation 



Isc( SS ) ((2ti-1)|,(2to -l)f) 



jnra 

sc(ss) 



will be helpful. 

Due to the properties: I SC {—X,Y) = —I SC {X,Y), 
I SC {X,-Y) = I SC (X,Y), I SS {-X,Y) = -I 3S (X,Y), 
I SS {X,Y) = I SS (Y,X), it is enough to define I™ 1 and 
J" s m only for n, m being positive nonzero integers (natural 



numbers), thus they can be presented as infinite matri- 
ces. In the thermodynamic limit, one has (see Appendix 
A for the general expression) J" c m 
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Note that there was no reason for presenting I, 



(10) 



sc 1 ss 



as matrices, other than the convenient visualization. 

We have found that the exact values of ip r , provided 
by |8]), (O, (fT0|) within the vortex core, are quite close to 
that of its asymptotic form found by Villain : 



Vv~5> R $ r (R) 



(11) 



R 



where $ r (R) is the polar coordinate of point r in the 
coordinate system with its origin at point R (the ref- 
erence angles are such that tp r — ip r r > ir if m r , r ' = 1, 
ip T — tpr' < —it if m r .r> = — 1, and |r/) r — ip r '\ < tt if 
m rr ' = 0), and 



9r = mi 2 + m 2 , 3 + m 3 4 + m 4j i 



(12) 



is the topological charge defined at site R of the dual 
lattice (see Fig. [1]). 



O 



R 



FIG. 1: Plaquette of sites 1,2,3,4 of the initial lattice adja- 
cent to site R of the dual lattice. 

Compare, for example, the field ip T given by Eqs. dHJ) 
and (|11[) for a vortex-antivortex pair with the minimal 
separation (sec Fig. [5]), shown in Tabs. Q] and [II] . 

It is worth mentioning that Eq. (|11[) can be derived 
from ©, using the asymptotic form of I sc and I sc , Eqs. 
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FIG. 2: Field ipnj) for a vortex-antivortex pair situated at 
sites (4.5.4.5) and (4.5,3.5) (the open and filled circles rep- 
resent the vortex and the antivortex, respectively) given by 
Eq. (fTTjl (Tab.|T]). The difference with the exact result follow- 
ing from (j8} (Tab. [TTJ) is insignificant within the resolution of 
the present picture. 



l9[) . (|A10[) . and integrating (instead of summing over 
R) along a properly chosen path L (or along several paths 
for the vortices with higher topological charges, see Sec- 
tion IA) connecting the vortex with its antivortex. We 
have verified that the form of the field ip T given by Eq. 
([5]) is independent of the particular form of this path L. 



TABLE I: Field ip^ (see Fig. [2) given by Eq. JTT]) 





i = 1 


i = 2 


i = 3 


i = 4 


3 


= 1 


0.1651 


0.1651 


0.1355 


0.0555 


J 


= 2 


0.2154 


0.2450 


0.2450 


0.1244 


J 


= 3 


0.2630 


0.3430 


0.4636 


0.4636 


J 


= 4 


0.2838 


0.3948 


0.6435 


tt/2 



B. Interaction between vortices in the Villain and 
2D XY models 

If i/j T is given by Eq. ([8]). the Hamiltonian ([6]) can be 
reduced to 

H = \ E (^r - fr'f + 9RqwV(K - R) , (13) 

<r,r') R,R' 

where topological charge q T is defined by Eq. (fT2)) . Now 
the vortex interaction energy is given by the second term 
in the Hamiltonian (1131) with 



V(R - R') = ^^£^. fc -( X - X ') cos ^( y - y ') 



2 k x a 



2 k v a 



(14) 

In the thermodynamic limit one can replace the sum 
over the first Brillouin zone in (|14j) with an integral, and 
then, since the difference between the Cartesian coordi- 
nates of the vortices centered on sites of the dual lat- 
tice is always an integer number of lattice spacing a: 
X — X' = na, Y — Y' = ma, following the same scheme of 
integration which was applied in Appendix A to obtain 
H}, El]) , one has: 



(-l)*(2n)! j2 (-!) J (2m)! 



V(na, ma) - g ^ _ ^ ^ . )m)] 



{-l) k {n-i)\ ^ {-l)\m- 3 )\ 



(n. — •) 



E 



(n-i-fc)! k\ ^ (m-j-l)\ 11 

k=Q y ' 1=0 y J ' 



F(i + k,j + l) 

(15) 



with F(p, q) given by Eq. (|A5|I . 

Then the energy of a vortex-antivortex pair, q-^ = +1 
and gR' = —1, which follows from (|13[) . is E pa \ r (x) = 
— V(x), where x = |R — R'| is the distance between the 
vortex and the antivortex. Comparing -E pa i r (x) that fol- 
lows from (fl~5)) with the asymptotic expression found by 
Villain 0: 



-Ep air (x) ~ 10.158 J + 2irJ\n{x/a) 



(16) 



TABLE II: Field (see Fig. [2} given by Eq. 



J = l 
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i = 4 
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n 7T 
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63^ 296 13 _ nn 3„. 
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0.7124 



j = 4 



118 _ 25 
5 2 

~ 0.0635 



71 



0.1460 



8-fTT 



2 

~ 0.4292 



tt/2 



see Fig.[3l we notice a fine agreement. The low number of 
points for small x/a is due to limited number of possibil- 
ities to situate the pair on a lattice, and the "oscillation" 
of data is the anisotropy effect for different orientations 
of vector R — R'. 

As it was already mentioned in the Introduction, the 
cores of vortices are characterized by large angles between 
the neighboring spins, so the harmonic approximation ([3]) 
cannot give the correct value of the energy of vortex cores 
in the model ([T]). Obviously, this can lead to different 
inter-vortex interaction energies in the Villain and 2D 
XY models. 

To estimate the energy of the vortex-antivortex inter- 
action in the 2D XY model, E%%? Y (x), we consider field 



< = $ r (R) - $ r (R') 



(17) 



5 




Eq. (15) 
Eq. (18) 
Eq. (16) 
Eq. (19) 



ln(x/a) 



FIG. 3: The energy of a pair of topological charges qn = +1, 
q R i = — 1 in the Villain and 2D XY models as a function 
of the separation x = |R — R'|. Open circles represent the 
exact result for the Villain model, following from (I15|) ; filled 
squares represent the numerical result for the 2D XY model 
(see HB])); solid and dashed lines are the asymptotic forms 
<Tl6t and JT9i 



(see Eq. (fTTj) ). which corresponds to the topological 
charges (/r = +1 and q^y = — 1. and assume that 



p2DXY i 



£ 

(r,r'> 



(x) = JY] [1-cos(^-^)] 



(18) 



performing the summation numerically over a system of 
sufficiently large size. We are aware that this assumption 
is not grounded, since ip r and ip r cannot be decoupled in 
the Hamiltonian ([T]), but it may be instructive. 

The quantity which is accessible for measurement in 
Monte Carlo simulations is the vortex-antivortex pair cre- 
ation energy in the 2D XY model, i.e. the energy of a 
vortex and its antivortex at the minimal separation a (see 
Fig. [J): E^* Y (a). The microcanonical Monte Carlo 
simulations showed that E^* Y (a) ~ 7. 3 J 0], while 
the canonical MC simulations gave 7. 55 J Q. Our es- 
timation which follows from (fT5|) is E^f Y (a) ~ 6. 6 J, 
in reasonable agreement with the mentioned computer 
experiments (the exact result for the Villain model is 
tt 2 J ~ 9.9J, see (fT5|) ). 

Comparing the result of (fT5|) to the vortex-antivortex 
interaction energy in the Villain model, see Fig. [31 we 
see that while at large separations 



III. INTERACTION BETWEEN VORTICES 
AND SPIN VACANCIES 

A. Hamiltonian of the Villain model with spin 
vacancies 

With the use of variables c r , taking values 1 and 
depending on whether site r is occupied with a spin 
or "empty" , respectively, the Hamiltonian of the Villain 
model with spin vacancies can be presented as (22j : 



H 



J 



8 r i — 27rm r iF /) c r c r 



(20) 



<r,r') 



Alternatively, it can be written via variables p r 
as 

H = H +£ft-#i(r) + PrPr>H 2 {v,v') 

r (r,r') 



(21) 



where Hq is the Hamiltonian of the Villain model without 
vacancies (Eq. (f2T)]) with all c r = 1), 



J 



Hi(r) = --J2(0r- 



? r+u - 2nm r , r+u ) 2 (22) 



with u = (±a,0), (0,±a) is the change in energy caused 
by the removal of the four bonds adjacent to the spinlcss 
site r, and 

H 2 (r,r') = ^(0r~ 9 r < - 2^m r ^) 2 (23) 

compensates the double removal of a common bond of 
two neighboring sites r and r' when there happen vacan- 
cies on neighboring sites. 

Applying ([5]), one can distinguish in the Hamiltonian 
(|20l) terms dependent on vortexless field ip r and/or vortex 
field ip r (marking them with indices tp and ip): 



(24) 



The first term, considered separately, describes a system 
of planar spins with angles tp r on a diluted lattice in 
the SWA, which was the subject of studies [l5| and (23| . 
for example. Here, we focus primarily on the last two 
terms that are connected to the presence of vortices in 
the system. 

Notice that taking tp r in the form of Eq. ([8]) does not 
lead to decoupling of ip r and ipr in the Villain model with 
spin vacancies (H ip ^ ^ 0). 



E 



2DXY 



(x) ~ 8.1J + 2irJln(x/a), 



(19) 



B. Hamiltonian of the Villain model with spin 
vacancies in the Fourier-transformed variables 



Ep®* Y (x) deviates considerably from the logarithmic 
form as the vortex and its antivortex approach each 
other. 



Fourier transformation of variables ip r , ip T and m r y 
allows to manipulate Hamiltonian (|20[) with much ease. 
The corresponding Fourier transforms ipu_, ipk and to£ 
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(a = x, y stands to distinguish two sets of Fourier trans- 
forms that correspond to "vertical/horizontal" orienta- 
tion of bond (r, r')) can be introduced via the following 
relations: 



0k 



^£ e -v, ^ = ^£*- ik V k! (25) 

r k 

4fE e ^- ^ = ^=J2e-^, (26) 



v r 

1\7 £— < K ' 



(27) 



a = x, y , 



where N is the number of sites in the lattice and the 
sums over r and k span the original lattice and the 1st 
Brillouin zone of the reciprocal lattice, respectively. Note 
that in (|2"T|) = (a, 0), u y = (0,a), so the property 
n^r,r' = —mr'.r is supposed to be used to obtain the 
Fourier transform of m r+VL r . 



J 



Then, for the field ijj T given by Eq. (|8|) one has the 
Fourier transform Q 

ml sin k£ + ml sin ^ 
V'k = -™ 2 , 2 \ ka 2 ■ 28 

Now, using f28|) . the condition ([7]) can be easily checked. 

After applying (|27|) and introducing the Fourier-trans- 
form of the topological charge q r : 



<7k 



-ikR 



<Zk, (29) 



Eq. ([12]) takes the form 



Then, it is quite straightforward to obtain (the reader 
is referred to Eqs. (|2~Tj) and (|2~i|) to understand the upper 
and bottom indices in the left sides of the equations) 



w , 4ttJ cos 
#i ( r ) = "TT E E ^ 



(k y +k')a . k y a ■ k'a 

. sm _|_ sui _|_ _ C os f*— sin -J- sin - j- _ i(k+k> 



k k' 



sin 2 ^£ + sin 2 k£ 



(31) 



t^j cos (^^ s i n ^ sin ^_ C0S (^Z±^> sin Ysin^ 



k k' 



/ • 2 fc~a i • 2 k v a\ ( . 2 fc~a . . 2 fc i a 

I sm =2 + sin -|-J Ism -|- + sm -|- 



i(k+k')r 



(32) 



^^^k'e * (k+k)r <V_ r . Ul e 1 ^ +^'-r.-u x e l " )sin¥sin-» 



N 



k k' 



5r'-r,u„ e 4 ^ + (5 r '-r.-u„ e* 2 



)• fc„a • k'a \ / . 2 i ■ 2 fc y a 

sm -|- sm -|- I sm -|- + sm -|- 



(33) 



Hf(r,v') = -^EE»' e " ,(k+k ' )r ( (^-.» 



2N 



k k' 



r — r,u„ 



J r' — r,-u„ 



jj . k v a . k' y a 

<Sr'-r.-u, e l 2 | sm -|- sm -§- 



3in Y ¥ ) / E sin2 V E sin2 ¥ ■ ( 34 ) 



a—x,y 



ot—x^y 



C. Attractive interaction between spin vortices and a spin vacancy 
Returning to variables ip r and in Eqs. (|3"Tj) , (j3"2")l . one has 

Hf^(r) = 7rj£q R { (<^ r+Ux - v? r _ u J/ sc (y- Y, x - X ) - (^r+u„ - ^r-u») 4c(z - X, y - Y) 

R 

+ (fr+u y +Pr-u y ~¥r+u x ~ <^r-u x ) ha ~ X, y - Y) } 



(35) 
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and 



Hf(r) 



k 2 JJ2Y1 «a«a'{^e(a; ~X,y- Y)I sc (x -X',y- Y') + I sc (y -Y,x — X)I sc (y — Y',x— X 1 ) 



R R' 

+ 2I ss (x-X,y-Y)I ss (x-X',y-Y')} 



(36) 



where I sc and I ss arc defined by Eqs. (|Alj) and (|A2|) . 
Analogous expressions for (|33p and can be obtained 
easily. 

In order to obtain the effective Hamiltonian describ- 
ing interaction between vacancies and topological charges 
only, one has to integrate out ip r in the partition function 



so that 



(37) 
(38) 



where Effa is the desired Hamiltonian. 



We have to restrict our consideration to the case of 
one spin vacancy at site r* to be able to use the results 
of Appendix [B] Then, using (|B2[) and (|3~Tj) , one has the 
effective Hamiltonian 



(tt-2 

k^CL • k'a - fc„ a 

-§- sin sin -|- 



k,k' 

- 2-^—^ sin sin ^ 
tt-2 2 2 

Finally, using ([3"fJ| and (j2"9")l . one can write: 

R R' 



z„a k v a • A:' a k' v a , • La fc^a • ^'a k'a 

-Tj± cos -|- sin -|- cos -| — h sin -|- cos -|- sin -|- cos -|- 
e -i(k+k>y ^ ^ s in 2 % 



ct—x,y 



(39) 



7 sc ( 2 /*-F,x*-X)7 sc ( 2/ *-y'. 



X')} I ss (x*-X. 

7T — I 



y* - y') 
y* - Y)i m 



x* - X 



,y*-Y')}.(40) 



While ([A3]) and (|A4|) provide the exact value of (gOJ) for 
a discrete lattice, it is instructive to find its asymptotic 
form: 



fl&(r*) = -(7r-l)J a 2 ^g R ^ 

R,R' 

r* -R)(r* -R') „ / 



R| 2 |r* - R' 



O 



\r* - R| 2 |r* 



R' 



(41) 



which follows from (|A9|) and (| A10|) . If, for example, one 
has a vortex of topological charge either + or — 1 and 
a spin vacancy, separated by distance x, (|41[) gives the 
energy of their interaction: 



E{x) = - J(tt - l)a z /x z + 0(l/x z ) 



(42) 



(compare it to the exact result following from Eq. (QHJ) 
shown in Fig. @|. 



D. A vortex pinned by the vacancy 

An analogue of the condition (JT]) for the field tp r , which 
would assure that = in the diluted Villain model 

(|2H|) . reads as 

c r (ip T — V'r+u — 27rm riI . +u ) c r + u = for all r . 

(43) 

Numerical studies of spin vortices in the presence of 
a spinless site [lcj . [Tl1 | suggest that it is energetically 
preferable for a vortex to be pinned (centered) on the 
vacancy. Thus, one can assume that 



±$ r (r*) 



(44) 



where < I ) r (r*) was defined after Eq. (fTTj) and r* is the 
coordinate of the vacancy, might satisfy ([4"3"f when the 
topological charge q = ±1 is on one of the four dual 
lattice sites R* adjacent to r* (see Fig. [5]). 
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-2 - 



i -3 



analogous to that of gSJ) leads to E*P XY ~ -3.21J, 
which agrees with —3. 178 J of the energy minimizing it- 
erative method and — 3.54J of the spin dynamics 
simulations [ill. 



IV. CONCLUSIONS 



Eq. (40) 
Eq. (42) 



0.5 1 



1.5 



2.5 3 

x/a 



3.5 



4.5 



FIG. 4: Interaction energy of a vortex of charge ±1 and a 
vacancy as a function of their separation x. Open squares 
represent the exact result (|40[) and the solid curve represents 
the asymptotic expression ()42[) . 



^ A r* ^ 

| R * 



p 7* T* 



FIG. 5: Representation field ip T of topological charge q = 
+1 situated at site R* of the dual lattice, which leads to 
its decoupling in the Hamiltonian of the pure Villain model 
(black and white arrows) and a model with a vacancy at site 
r* (grey arrows). 



Then, the vortex-on-vacancy pinning energy, i.e. the 
energy of the vortex centered on r* minus that of the 
vortex centered on R.*, can be estimated by numerical 
summation over a lattice of sufficiently large size: 



-Epin = — ^2 [ (y> r - 1pr> ~ 27TTO r 
<r,r') 

- {ip' r - i'' r , - 2nm r yY 



(45) 



where ip' r = ±<£> r (R*), which gives E ] 



pin 



-5.22J. 



The corresponding energy that follows from (|4TH is 
Spin = -(3tt - 4) J ~ -5.42J. The difference from the 

result of (|45p is not surprising, if one notices that ip r only 
approximately fulfills (|4"3")l for almost all the lattice sites. 

It is worth mentioning that using the exchange poten- 
tial of the 2D XY model, Eq. ([T]), numerical summation 



The exact and asymptotic expressions for the interac- 
tion energy of topological charges and a spinless site, Eqs. 
(|4U1) and (JUJ), found for the Villain model on a square 
lattice, definitively confirm the attractive character of the 
interaction. This agrees with the results of the spin dy- 
namics simulations for the 2D XY model and the 
energy minimizing iterative method for the easy-plane 
Heisenberg model [Io| . 

However, we showed that this interaction in the 2D 
XY model can differ from Eq. (|4T)]) . that corresponds 
to the Villain model, considerably at small separations 
due to different energies of the vortex "cores" (regions 
with strong disorientation of spins). In particular, the 
exact value of the vortex-on-vacancy pinning energy in 
the Villain model, E pin = -(3tt - 4) J ~ -5.42J, differs 
significantly from that found in [ll[ and [l(| (—3.54 J and 
—3. 178 J, respectively). 

Moreover, we showed that the mentioned difference of 
the vortex cores' energies in the two models leads to a 
deviation of the vortex-antivortex interaction energy in 
the 2D XY model from a logarithmic law at small sep- 
arations, while the corresponding energy of the Villain 
model retains logarithmic dependence on separation x 
(if we neglect slight anisotropy effects) up to the small- 
est possible distance on a lattice which is of one lattice 
spacing, x = a. 

We have estimated the vortex-antivortex pair creation 
energy for the 2D XY model as E^f Y (a) ~ 6. 6 J (in 
contrast to ir 2 J ~ 9.9J of the Villain model), which is 
in reasonable agreement with the results of the recent 
Monte Carlo simulations 0, (7.55 J and 7.3 J, respec- 
tively) . 
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Appendix A: Functions I SC (X, Y) and I BB (X,Y) 

In this appendix we study the functions 
1 



I SC (X,Y) 



N 



E 



sin %2 C os %2 

sin k„X cos k„Y 



sin sin _a£ 

2 fc Q a 



(Al) 



sin/c^A sin/c,.Y , 



which enter many important expressions concerning the 
behavior of topological charges, and the sums over k in 
(|Alj) . (|A2|) span the first Brillouin zone. 



For X = (2n — l)a, Y = (2m - l)a (n, m = 1,2,3,.. .) 
()Alj) and (|A2|) can be calculated exactly, replacing the 
sums with integrals in the thermodynamic limit. The 
(A2) integration gives 

I 



J sc ((2n- l)o,(2m- l)a) 



(-lp'(2m - 1)! (~l) k (n - * - 1)! 



V- 1 (-l) j (2n- 1)! ^ 

^ (2(n-i-l))!(2i + l)! ^ (2(m - j) - l)!(2j)! ^ (n - i - * - l)!fc! 



•Til E 



m-j 

E 

i=0 



(-l)'(m-j)! 
(m-j - Z)! Z! 



(A3) 



, //„ ^ / n w v- 1 (-1)*(2»-1)! ^ f-l) J *(2m-l)! 
7„((2n- l)a,(2m-l)a) = > — — ; v sst , ' > 



E 



(-l) fe (n- i- 1)! 



Q (2(n-i-l))!(2i + l)l ^ (2(m-j-l))!(2j + l)! ^ (n-i-*-l)!fc! 

x m E 1 m + * + i,i + 1 + 1) (A4) 



(m-j -I- 1)! Z! 



with 



= «i (2(p + u)-l)II !)-!)!! 1 P ^ (-!)«+> + g - l)\(2u 1)11 

[P,q) ' (2(p + «))H (2(g -«-!))!! 2^ Q (p + g - « - 1)! (ulf 



<7T 



+ 1) / ( 2u - l)!! " 1 (u-to-l)! 1 



7T ^-^ (p + g — u — 1)! it! 



E 



u! ^ (2(u-w) - 1)!! u 
I 



(A5) 



These results were obtained by expressing smk x X, 
sin k y Y, cos k y Y as polynomials P (sin %p , cos , and 
then applying the standard tables of integrals [25(. We 
used the notations: (2n)U = Y\7=i (2"- — 1)!! = 
nr=i(2i-l); whenn = 0: (2n)H = 1, (2n-l)!! = 1. The 
sums of no meaning, like _ with m < n, that may 
be encountered in (|A5[) for some values of p, q, should be 
interpreted as equal to zero. 

It is instructive to find an asymptotic form for (|A1[) and 
(|A2[) . It turns out that simple analytic expressions can 
be obtained, assuming that at least one of the arguments 
X, Y is large. Using the integral [24j 



one can show that 
I sc (A^oo,Y) 



pn/a 

-J dk y e~ x ° sm — cos k y Y 

■ i I ■ K y' 
x sinn sin — 

V 2 , 

a r^ a 

- \ dk y e~ xkv cos k y Y (A7) 
n Jo 



k v a\ 

cot 



kyCl 

~2~ 



and 



2\a\ 



-\a\ 



(A6) 



f'-K I a 

I SC (X,Y ^ oo) = - dk x e~ Y ^ sin ^ cosk x X 

I" JO 

dk y e~ Ykx sink x X. (A8) 



/ . k x a\ 




r— J 






I" Jo 
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So, 



a X 
I SC (X,Y) = — 



7T X 2 + Y 2 ' 



(A9) 



when at least one of the arguments X, Y is sufficiently 
large. 

In a similar way one can show that 



I SS (X,Y) 



a 2 XY 



(A10) 



7T {X 2 +Y 2 Y 

if at least one of the arguments X, Y is sufficiently large. 



Appendix B: Hamiltonian describing topological 
charges in a system with a spin vacancy 

The aim of the present appendix is to show how the 
"vortexless" degrees of freedom ip T can be integrated out 
in the partition function (|37p when only one spin vacancy 
at site r* is considered, so that (see Subsection IIII A|) 



0, r = r*, 
r^r*; 



or p r 



1, r = r*, 
r ^ r*. 



(Bl) 

As we will show below, the partition function can be 
presented in this case in the form (|38[) with 



ME 



<7k<?-k 



4/3 2 Jy — 7k 



7k 



V-UTj-k^ -j(k+k')r* 



7T 1 \ - 



k,k' 



v) V-*V-u e -iQt+X)r- ] (B2) 

7k' J 



where 7k = 2 (sin 2 %2 + s i n 2 _|_^ ; 

3k,k' = (7k+k' - 7k - 7k')/7k, (B3) 

and 



»7k 



* u 

(B4) 



1. The partition function of the Villain model on 
the lattice with a spin vacancy 



Let us denote 



so 



Z. 4 , = Tr v e- piH ^ +H '<'-*\ (B5) 

(B6) 



Using Fourier transformation (|25[) , one can rewrite the 
terms that depend on ip v in the Hamiltonian (|24[) as (see 
El): 



k,k' 



(B7) 

where gk,k' was defined in (|B3[) and the sums are over the 
1st Brillouin zone. Correspondingly, the mixed tpijj-term 
in Eq. reads: 



— ikr* 



x ]T e-* ku (^ r , +u - ^ - 2^m r . +u , r . ) j . (B8) 

u / 

Using the Taylor series expansion, Z^ can be written 



as: 



Zih = Tr e ~' 9J ^-' k7k¥ ' k ' p ~ k+ ^ kJ,kVk 



/ oc 1 

1 + E^ 

\ 71 = 1 



v k2 ),-(Vk 2 „_ 1 ,y 1 c 2 „) (B9) 



where 77k was defined in (|B4[) and 

_ {-PJ) 



'(Vkj ,fk 2 )-(*'k 2 „_ 1 >¥>k 2rl ) 



^E- E 



ki,k2 kl2 n — l,k2n 

x e-'' kl+ ^ k2 »» r ' ffkl;k2 ... 9k! „_ l!k ^ kl .^ k2 , (BIO) 
Now, introducing the notations: 



Z„ = r Yr v e^ l3J ^ k ' lkipk ' p ^ k+ ^ k '' ,k ' Pk (Bll) 



and 



- 1 Tr v (e 



~PJ Ek 7k¥>kV?-k + X] k Ik^k 



(B12) 



B9D becomes 



(00 ^ 
1 + E a ( 
n— 1 



(B13) 



2. Calculation of Z, and ((^k^ • ■ • V?k 2 „)* 



Z» and (yki ••• ( Pk 2 „)* are t ne nrs t quantities to be 
calculated. 

Since tpw (for k ^ 0) is a complex quantity: ip^ = 
tp^ + if l--, Tr v should be understood as: 

/oo poo 
dp£ / dy£ , (B14) 

where B l / 2 stands for a half of the 1st Brillouin zone 
excluding k = (y> k and </? k in the other half are not 
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independent, due to the relations: ip c _^ = <p k and y?f_ k = 
— ffj. It was possible to extend the bounds of integration 
to infinity in (|B14[) and omit writing the integral over 
</?0j since the functions that stand after the trace in our 
calculations are always rapidly decaying when (3 J — > oo 
and independent from <po- 

Then, it is straightforward to obtain 



n 



7T \ 1_ V ''k'i-k 



2/3J 7r 



(B15) 



Using (|B15|) . it is easy to show that 



Z* , (B16) 



dr] kl drjk 
where 

d d . d d d . d 

Noting that = 2<5 kjk < (#k,k' is Kronecker delta), one 
arrives at 

n 1 

;=0 ^ ^ ' pairs 2n->Z 



n 



2n-2l 



3k,„,-k Jt , T-r >?-kp„ 

7k,„ 11 7k„ 



, (B18) 



where the sum spans all the possible ways of 

pairs 2n— 

selecting I indistinguishable unordered pairs (z u , j u ), u = 
1, . . . , I out of 2n indexes 1, . . . , 2n. (It is easy to see that 



E i 



(2n)! 



pairs 2n— >■/ 



i!(2!) J (2n-20! 



3. Calculation of ( J 



(<Pk 1 >V>k 2 )>-'-(¥>k 2 „_ 1 ,<f>k 2 „) 



According to (|B10|) and (|B12|) . 



(Vkl <¥ 3 k 2 ),---(¥>k 2ll _ 1 ,¥>k 2 „) 



(~/^) r 
N n 



E 



kl,...k2n 

xe -i(k 1 + ...+k 2 „)r* Ski ^ . . . flkan _ i)kaB . . . <^ k2n ) t . 

(B19) 

At this stage, it is convenient to introduce the notions: 



II - TV' 1 E 



5ki,-k 2 5k 2 ,-k 3 ' ' ' ffki-i.-kiffki.-ki 



ki k; 



(B20) 



ki,...,k i+ i 
i(ki+k i + i)r* ^ki^kj + i 

7k l+1 



x e 



(B21) 



Then, insertion of (|B18[) into (|B19[) leads to a polyno- 
mial form with respect to Ii and I* (i = 1, . . . , oo): 



(Vkj >¥ , k 2 )>---(¥ ! k2 n _i >¥>k 2 „) 



= (-i) n E 



1=0 



(2/3J)' 



Z fn-l[(n-l)/j]\ / I \ 

■ nz n e p ^-' 

^=lAi=0/ y=i A*=0 / \i=l 



(B22) 



where [a] means the nearest integer not exceeding a, 

1, £ = 



<*(*) = { 



0, x ^ 



and 



iV Ai,...,A, 



"11 Aj!(i!) Ai 11 A*!(z!)^ 



(B23) 

is the combinatorial "weight" given by the number of 
ways of selecting Ai unordered elements, A2 unordered 
groups of two unordered elements, A; unordered 
groups of I unordered elements, A* unordered elements, 
Aj unordered groups of two unordered elements, and 
A* _j unordered groups of n — I unordered elements out 
of n distinct elements, which is 

n!/(Ai!A2!"-A,] A*^! ■ ■ • A*_ ; ! (1!) A ^(2!) A = ■ ■ ■ (Z!)^ 
(l!) A i(2!)*S...((n- 00^-0 , 

times the number of distinct ways of connecting 4 distinct 
elements belonging to 2 distinct groups, each consisting 
of two elements, with 2 indistinguishable links in such a 
manner that the two elements of one group are connected 
to the elements belonging to another group, raised to the 
power A2, times the product over i = 1,1 of the number of 
distinct ways of connecting 2i distinct elements belonging 
to i distinct groups, each consisting of two elements, with 
i indistinguishable links in such a manner that the two 
elements of each group are connected to the elements 
belonging to two another groups, raised to the power Ai, 
i.e. 



n^-i)!]' 



(B24) 



times the number of distinct ways of connecting 4 distinct 
elements belonging to 2 distinct groups, each consisting 
of two elements, with 1 link in such a manner that one 
of the two elements of one group is connected to one 
of the two elements belonging to another group, raised 
to the power A2, times the product over j = l,n—l 
of the number of distinct ways of connecting 2j distinct 
elements belonging to j distinct groups, each consisting 
of two elements, with j — 1 indistinguishable links in such 
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a manner that one of the two elements of any group is 
connected to one of the two elements of another group 
and the second element is either connected to one of the 
two elements of a different group or not connected, raised 
to the power A*, i.e. 

n—l 



II y 

0=1 



- 1 -.'H Aj (B25) 
Inserting (|B23|) into (|B22|) . one has 

n 1 
(-!)"»! £ ( 2 ^J)n-i 



(Vkj ,Vk 2 ). — (Vk an _l .¥>k2n) 



1 = 



I [l/i] . , T x A* n-/ [(n-0/i] - / r *\ AI 

xIIEi(l) II E ^ | 

i=U,=0 v 7 j = l A*=0 1 x 



1 



n—l 



x S ( J2 iX i - 1 ) S \ EA* - (»-0 ) . (B26) 

xi=l / 



and then, inserting (|B26I) in (|B13|) . one can notice that 
the infinite series in (|B13I) can be rearranged as it is 
shown below: 

z, = z.n(i + (-D'|4((-ir|) 2 
4 MD^-) n ( i+ <- i » ) w 
+ K<- i ) , l;) i+ 5!( ( - 1)J w) 3+ ^)' 

or 

^ = ^exp^^-l)^/^ 

x-P^B- 1 )^- (B27) 

4. Calculation of /; and I* 
Eqs. (|B20[) and (|B21|) can be written as 

^ = ^EW k >- k ) (B28) 

k 

and 

7 * ^E^^ 14 ')^ e^ k+k >* (B29) 
(i > 1) with 
'>> k ')-^ E 



ffk,-k!3ki,-k 2 ' • • ffki-i.-kiffk^k' 



ki,...,ki 



for i > 1 and Io(k, k') = 3k, k'- One can notice the obvi- 
ous recurrent relation 



/ i+1 (k ) k') = i^/ i (k,-k , ) 5 k.,. 



(B31) 



In the thermodynamic limit, one can replace the 
sum ^ k over the 1st Brillouin zone by the integrals 

J^Y f-i°/a dk * J-i°/a dk v> and tncn > noticing that 



a 2 r /a r /a sin 4 %^ 

- / d H *% in2 ^ +sm ^ 



7T 



1 

7T 



and 



a 2 r/ a sin 2 ±f cos 2 ^ _ a 2 

^ dk X ^ dky-^-——^ - - 

r n/a sin 2 fc^a sin 2 1 j 



/0 Jo 



2/ ^ 2 fe | S +sin 2 ^ 2 7T ' 



one can show that 

^E^-k*,k» = ( 1 - 7 ) - 9k '- k " 



-- (Sk,-k" + £fk,k") + ( o ~ _ ) 

7T \ Z 7T 



3k,k" 



-- (gu-k" +5k,k") + ( - — j 

7T \2 7T, 

and 

E 7k5k,k' = -7k' ■ 

k 

Then, it is easy to see that: 

/,(k, k') = A i 5 ki (_ 1) i k / + 5, (3k,-k' + 3k, k') + CiTk' 

with coefficients A4, Bi and Cj, obeying the recurrent 
relations: 



= ( 1 - - Mi , 



B i+1 = --Ai + (l--)Bi, 



Ci+1 



1 1 



(B30) 



, 2 ^ 

and A = 1, B = 0, C = 0. Thus, 



At = I 1 - - 

7T 



(A; + 25, ) - Ci 
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Finally, one can obtain expressions for and I* and 



check that 

oo 

V(-i)' 1 /* = - — 

£f 4(tt-2) 

X T? Z^(-9k -k' +5k,k') e ^ ; 

k,k' ,K 

+ l4E(^-^ +Sk,kO !L ^ e - 4(k+k> *-(B32) 
4 JV 7k' 
k,k' 

In conclusion, using ([B"6)l. (|B15|) , (|B27|) . and (|B32|) . we 
obtain (|B2l . 



[1] N. D. Mermin, Rev. Mod. Phys. 51, 591 (1979) 

[2] P. M. Chaikin, T. C. Lubenskii, Principles of condensed 

matter physics (Cambridge University Press, Cambridge, 

1995) 

[3] V. L. Berezinskii, Sov. Phys. J.E.T.P. 34 610 (1972) 
[4] J. M. Kosterlitz, D. J. Thouless, J. Phys. C 6, 1181 (1973) 
[5] J. M. Kosterlitz, J. Phys. C 7, 1046 (1974) 
[6] J. Villain, J. Phys. (Paris) 36, 581 (1975) 
[7] S. B. Ota, S. Ota, Phys. Lett. A 206 133 (1995) 
[8] R. Gupta and C. F. Baillie, Phys. Rev. B 45 2883 (1992) 
[9] L. A. S. M61, A. R. Pereira, A. S. T. Pires, Phys. Rev. B 
66, 052415 (2002) 
[10] G. M. Wysin, Phys. Rev. B 68 184411 (2003) 
[11] A. R. Pereira, L. A. S. M61, S. A. Leonel, P. Z. Coura, 

B. V. Costa, Phys. Rev. B 68, 132409 (2003) 
[12] S. A. Leonel, P. Z. Coura, A. R. Pereira, L. A. S. M61, 

B. V. Costa, Phys. Rev. B 67, 104426 (2003) 
[13] G. M. Wysin, A. R. Pereira, I. A. Marques, S. A. Leonel, 

P. Z. Coura, Phys. Rev. B 72, 094418 (2005) 
[14] G. M. Wysin, Phys. Rev. B 71, 094423 (2005) 
[15] B. Berche, A. I. Farinas-Sanchez, Yu. Holovatch, R. Pare- 
des V., Eur. Phys. J. B 36, 91 (2003) 



[16] H. Chamati, S. Romano, Phys. Rev. B 73, 184424 (2006) 
[17] H. Chamati, S. Romano, Phys. Rev. B 75, 184413 (2007) 
[18] F. Wegner, Z. Phys. 206, 465 (1967) 
[19] J. Tobochnik and G. V. Chester, Phys. Rev. B 20, 3761 
(1979) 

[20] S. T. Bramwell, J.-Y. Fortin, P. C. W. Holdsworth, S. 
Peysson, J.-F. Pint on, B. Portelli, M. Sellitto, Phys. Rev. 
E 63, 041106 (2001) 

[21] J. V. Jose, L. P. Kadanoff, S. Kirkpatrick, D. A. Nelson, 
Phys. Rev. B 16, 1217 (1977) 

[22] O. Kapikranian, B. Berche, Yu. Holovatch, Phys. Lett. A 
372, 5716 (2008) 

[23] O. Kapikranian, B. Berche, and Yu. Holovatch, Eur. 
Phys. J. B 56 93 (2007) 

[24] R. J. Tallarida, Pocket book of integrals and mathematical 
formulas (3rd ed., Chapman & Hall/CRC, Boca Raton, 
London, New York, Washington, D.C., 1999) 

[25] A. P. Prudnikov, Yu. A. Brychkov, O. I. Marichev, Inte- 
grals and series: Volume 1: Elementary functions (CRC 
Press, 1986) 



